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Let Sfe = -^{vik,. . . ,«jvfc)^, k = l,...,K, where {vik,i,k = 1, . . .} 
are independent and identically distributed random variables with 
£^11 = and Evfi = 1. Let Sk = (si, . . . ,Sfc_i, Sfc+i, . . .,sk), Pfc = 
diag(pi, . . . , pk~i,Pk+i,- ■ ■ ,Pk) and f3k = Pks]\ {SkPkSl + cr^I)"^Sfe, 
where pfc > and the l3k is referred to as the signal-to-interference ra- 
tio (SIR) of user k with linear minimum mean-square error (LMMSE) 
detection in wireless communications. The joint distribution of the 
SIRs for a finite number of users and the empirical distribution of 
all users' SIRs are both investigated in this paper when K and N 
tend to infinity with the limit of their ratio being positive constant. 
Moreover, the sum of the SIRs of all users, after subtracting a proper 
value, is shown to have a Gaussian limit. 

1. Introduction. Consider a symbol synchronous direct sequence code 
division multiple access (DS-CDMA) system with K users. The discrete- 
time model for the received signal y in a symbol interval is 

K 



where the is the symbol transmitted by user fc, sj. £ R is the signature 
sequence of user k and w G i?^ is the noise vector with mean zero and covari- 
ance matrix a^I. We also assume that the symbol vector x = (xi, . . . ,xk) 
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has a covariance matrix P where P = diag(pi, . . . with pk being the 
received power of user k, that is, Exj, = pk and that the symbol vector is 
uncorrelated with the noise (more details can be found in [13]). 

The engineering goal is to demodulate the transmitted Xk for each user. 
Assume that the receiver has already acquired the knowledge of the signature 
sequences. For user k, the linear minimum mean-square error (LMMSE) 
receiver generates an output in a form a^y where is chosen to minimize 
the mean-squared error 

(1.2) E\xk-sily\^. 

The relevant performance measure is the signal-to-interference ratio (SIR) 
of the estimate (see [13]), which is defined by 

(1.3) (3k=Pk4iSkPkSl + a^I)-^Sk, k = l,...,K, 

where Sk and P^ are obtained from S = (si, . . . ,sx) and P by deleting the 
A:th column, respectively. 

It is difficult to obtain clear engineering insights from (1.3) since it is 
dependent on the signature sequences. However, if signature sequences are 
modeled as being random, one may further proceed with the analysis using 
random matrix theory when the number of users K and the processing gain 
N approach infinity, that is, suppose 

V iV 

k = 1, . . . ,K, where {vik,i, k = 1, . . .} are independent and identically dis- 
tributed (i.i.d.) random variables. Rigorously speaking, if Vik are random 
variables, then (1.2) should be viewed as a conditional expectation and at 
this time it is also necessary to assume that the signature sequences are 
independent of transmitted symbol and noise. 

Indeed, considerable progress has been made in this area. For example, 
Tse and Hanly in [11] derived the asymptotic SIR under MMSE, a decorre- 
lator receiver and a match filter receiver and fluctuations of SIR have sub- 
sequently been considered in [10]. Some related results can be found in [13]. 
Also see [12], and references therein and see the review paper [1] concerning 
random matrix theory as well. 

However, there are still many open problems in this area. For example, 
Tse and Zeitouni in [10] asked: What is the empirical distribution of the SIR 
levels of the users across the system? Is this empirical distribution suitable 
for characterizing the asymptotic distribution of the SIR for a particular 
user? Is there any type of "weak asymptotic independence" among users? 
Also, the asymptotic distribution of the sum of all users' SIRs under MMSE 
has remained unsolved, which has a close connection with another important 
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performance measure, sum mutual information or spectral efficiency (suit- 
able scaling) (for more information concerning the sum mutual information 
or the spectral efficiency, see [9] and [14]). 

In this paper we will answer the above questions. In other words, we will 
derive the joint asymptotic distribution of the SIRs for different users and 
the limiting empirical distribution of the SIRs of the users across the system. 
The sum of the SIRs for all users, after subtracting a proper value, is also 
shown to have a Gaussian limit, which gives the asymptotic distribution for 
sum mutual information under MMSE. 

Before stating our main results, we will introduce some notation. Write 
Btv = SPS"^, whose empirical spectral distribution (ESD) is denoted by 
pBN The ESD of power matrix P is denoted by Hn- Let cn = K/N. 
F'^'^{x) and H will denote the weak limits of the distribution functions 
p^N ff^ respectively, as N, K ^ oo if the limits exist. Define b = J{x -\- 
a^)-^ dF""'^ (x) and In = /(x cj^)"^ (iF^^'^^(x), where F'''^'^'^ (x) = 
F'''^{x)U=,^^H=H^. 

Theorem 1.1. Suppose that: 

(a) {vij,i,j = 1, . . .} are i.i.d. with Evu = 0, Evfi = 1 and Evfi < oo. 

(b) Hfyf converges weakly to some distribution function H and the ele- 
ments of P are bounded by some constant. 

(c) K/N ^OO as N ^oo. 

Then, for any finite integer m 

(1.4) ( ViV(/3i - pibN), . . . , VN{f3m - PmbN)) ^ iV(0, C) 

with covariance matrix 

(1.5) 

X dmg{pl,...,pl,). 

Remark 1.1. Theorem 1.1 indicates that the asymptotic independence 
of the SIR among users holds, as conjectured by Tse and Zeitouni in [10]. 
This theorem also includes Theorem 4.5 of [10] as a special result. Actually 
Tse and Zeitouni in [10] only derived the asymptotic distribution for a single 
SIR under the conditions that pi = ■ ■ ■ = pk and vu is symmetric. 

Theorem 1.2. In addition to the assumptions (b) and (c) of Theorem 
1.1, we suppose: {a') {vij,i,j = 1, . . .} are i.i.d. with Evu = and Evfi < oo. 
Then 



(1.6) 



GNix)^G{x), 
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where 

1 ^ 

(1.7) GNix) = -Y.IiPk<x); 

k=i 

i.p. denotes the convergence in probability. Moreover, the Stieltjes transform 
ofG{x) is j\bx- z)-^dH{x). 

Remark 1.2. Theorem 1.2 characterizes the empirical distribution func- 
tion of the SIRs for different users, and, simultaneously, it reveals that the 
asymptotic empirical distribution of the SIRs for a whole system is different 
from the asymptotic distribution of the SIR for a particular user, which is 
normally distributed, as shown in Theorem 1.1. For example, consider a spe- 
cial case pi = • • • = pk = V'l then one can easily obtain G{x) = I[pb <x< oo) . 

Remark 1.3. Indeed, the convergence mode in Theorem 1.2 can be 
strengthened to converge with probability 1 according to Theorem 7.1 in 
[4]. In that paper, a more flexible model is employed; they show that the 
corresponding SIR converges with probability 1 and also provide uniform 
convergence of the SIRs for all users. It is interesting to consider how to 
derive the asymptotic distribution of the SIRs under their model. 

Theorem 1.3. In addition to assumptions (b) and (c) of Theorem 
1.1, suppose that {a") {vij,i,j = 1,...} are i.i.d. with Evu = OjEvfi = 
1 and Evf^ = 3. (d) J x{l + xb)-^dH{x) = JxdH{x)J{l + xb)~^dH{x) 
and Jx^dH{x)iJ{l + xb)-^dH{x)f + ] x^{l + xby^dH{x) = 2jx'^{l + 
xb)~'^ dH{x)f{l + xb)~'^ dH{x). Then we have 

K 

(1.8) Y.i(^k-bNPk)^N{i^,p), 



k=l 



where 



dH{x) 
{1 + xby 



, , , dF^^^ix) f x^ , , , 
1.9 =2c/- TT^dHix) 



' arg 1-c/- )-^dH{t)]dx 



and 



(1.10) 



27r7 (a;-Fcj2)2 J {l + tm{x)) 

dH{x) 



[i + xby 

1 f f {d/ dzi)m{zi){d/ dz2)m{z2) 



27r2 y J {z\ + a'^){z2 + a'^){rn{zi) — m{z2))'' 



■ dzi dz2, 
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where the contours for zi and Z2 are nonoverlapping and closed and are taken 
in the positive direction in the complex plane, both enclosing the support 
of F'^'^{x). Here m{z) represents the Stieltjes transform of F^'^(x) and 
'^m{x) = lim^^a; Qm{z) . 

Remark 1.4. Assumption (d) is satisfied when pi = ■ ■ ■ =pk = P and 
in this case the formulas (1.9) and (1.10) can be simphfied as 

' -^cfjm-J^dHix) 



(1.11) 



(l+p)'^ J {x + a'^/pYJ {l + xh) 
1 1 



+ 



4(a(c) + ctVp) 4(6(c) + a2^) 
1 dx 



and 

(1.12) P={l+P) 



27r Ja{c) {x + cr^/p) ^4c - (x - 1 - c)2 
2c 



((cj2/p + C-l)2+4cj2/p)2' 

where a(c) = (1 — a/c)^ and 6(c) = (1 -|- \/c)^ and the expression of F'^ is 
referred to [6]. 

From Theorem 1.3, we can obtain the fohowing corollary concerning sum 
mutual information under MMSE. 



Corollary 1.1. Under the assumptions of Theorem 1.3, 

K 

(1.13) E(log(l + Pk) - log(l + hNPk)) ^ N{iJii,pi), 

k=l 

where 

As is seen from Corollary 1.1, the sum mutual information normalized by 
N (which is the spectral efficiency which is relevant in wireless communica- 
tions) converges (see also [14]). So one can guess that the small fluctuation, 
when expanded by a factor of N, appears to be Gaussian. However, it is not 
an easy task to prove it. 

The organization of the paper is as follows. Section 2 establishes Theorem 
1.1. The proof of Theorem 1.2 is provided in Section 3. The proof of Theorem 
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1.3 and Remark 1.4 is included in Section 4 and the proof of Corollary 1.1 is 
contained in the last section. Throughout this paper, M may denote different 
constants on different occasions and || • || denotes the spectral norm of a 
matrix or the Euclidean norm of a vector. Also, set A = SPS"'" + cr^I, A^ = 
SfcPfcS|' + (T^I, k = 1, . . . , K, to simplify notation. 

2. Proof of Theorem 1.1. Before beginning with the proof, we first state 
a lemma. 



Lemma 2.1. Let ajj be the jth diagonal element of . Under the 
assumptions of Theorem 1.1, 



(2.1) 



1 



N 



hm — a,,: 



2 ^2 



Proof. From the well-known matrix inverse formula, we have 

1 



(2.2) 



where 



an 



(SiPiSf + C72l)n - sf Sf (SiSf + a2l)-iSis 



{^/P2Vj2,■■■,^/PKVjKy 



1 = (S2, . . . ,SArJ 



J = l 



, . . . , X » , 



A^, 



and (SiPiSj' + o-^I)!! is defined in (2.3). 

Applying the Helly-Bray theorem one can find 

K 



(2.3) (SiPiSf + c72l)n = ^f^Pfc<^ + a2^c / xdH{x)+a\ 

k=2 ■' 



where we also use the fact that 



E 



K 



k=2 



A^2 



k=2 



as A ^ oo. 

It is observed that 



Sj = diag(VP2, • • • , ^/VK)^j, Sf = diag(Y^, . . . , v^)Sf 

with Sj = ^j^{vj2, ■ ■ ■ :VjKy and Sf = (§2, . . . , stv)- This, together with 
Lemma 2.7 in [2], implies 

1 2 



E 



(2.4) 



s|S|(SiS| +cj^I)-'SiSi - — trPfSl (SiPiS{ +a^I)-^S 



< 



MEvi 



11 



N 



^||PiS[(SiPiSf + (t2i)-1SiP 



A^ 
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as — > oo, where we use 

||PiSf(SiPiSf + ct2i)-1SiPi|| 

<M||Sf(SiPiSf + (t2i)-1SiPi|| <M. 
Let s,- = {V2j,. . .,VNjf, = (SiPiSj + anr\ B^/ = (Si,Pi,S?;.+ a^)-^ 

with Sij and Pij obtained, respectively, from the matrix Si and Pi by re- 
moving the jth column. We then have 

4trp2sf(SiPiSf + a2l)-iSi 



N 



(2.5) 



i=2 

K 



K 



Pj 



i=2 



N 



i=2 



Hence it follows that 

K 



E 



1 K 

1 ^ Pj 



J_ 1 _ 

K 



1 



< M— y E 



1 



K 



(2.6) 



+ M—E 

N 



— trB7^ -6 

N ^ 



cK \ 



trB{ 



-1 



as ^ oo. 



In the last step, we also use ^trB;^ ^ h and the uniform integrability 
of N^^tiB^^. 

From (2.3)-(2.6) we have 

(2.7) sf Sf (SiSf + (j2l)-iSiSi j xdH{x) - c 
Thus combining (2.3) and (2.7) one can get 

(2.8) a..^(a^ + cl^-^V = h. 



X dH{x) 
l + xb 



l + xh 
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Since an is bounded by 1/cr , an is then uniformly integrable and so 
E 



1 ^ 



1 ^ 

as — > oo. Thus the proof of the above lemma is complete. □ 

We proceed with the proof of Theorem 1.1. Let Sik^ and Sifc^^j be the 
matrices obtained from the matrix Si by removing the fcith column, the 
/cith and A;2th columns, respectively and let Sim be the matrix obtained 
from Si by deleting the first m — 1 columns. The matrices Pifci , Pifcifco and 
Pim are defined similarly. Define A]"^^^ = (SifcjPifci Sf^^ +(T^I)"\ A.]^kik2 = 
{Sik.k^Pik.k^Sjk.k, + ^^I)"^ and A-^ = (Si^Pi^Sf^ + aH)~\ Further- 
more, all analogues such as A]'^"'^^^^^ required in the following derivation 
have similar meanings. 

Write 

pisf A^^si 

— PlSi^i^Si- 2^ — ^ 

(2.9) =p,s.A,,s,-E^^— j^^j— 

fc^2A:i^l (1 +PfclSriA^;,\ SfcJ(l +Pfc,S^^A-;-\;^^SfcJ 

^ Pl^l^lfnPk^^kAAlfnPk2^k2 ' ' • Pfc„,_i S^j^^^ A'^Si 

fc2^fcl,..., 

where the subscripts fci, . . . , /cm-i are larger than 1. 
For any j (i, j = 1, . . . , m) we have 

Ei^Kk^jf = EtrA-lsjsjA-ls,sJ = tr A^^ 



which implies 



Ar2 
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and then 

(2.10) VNsjA-lsk,slA-lsk, ^ 0. 
Hence, from (2.9) and (2.10) we have 

(2.11) \/iV(pisf A^^si - ViVpisf A^^si) ^ 0. 
Similarly, for any k= 1, . . . , m one can find 

^{pkslA^'s, - VNpkslA~lsk) ^ 0. 

It thus suffices to consider the asymptotic distribution for the linear com- 
bination of y/N (pfcS^A^^Sfe —pkbN), k = 1, . . . ,m. To this end, it can be seen 
that {pkS^A^^Sk,k = l,...,m} are independent when the matrix A]~^ is 
given and hence it suffices to consider only one of {pkS^A^^Sk, k = 1,. . . , m} 
when the matrix A^^J is given. 

By Lemma 2.1 and the Jensen inequality, it is easy to verify that 

N 

lim a^,- / a?- < 1, 
where A]~^ = (ojj). Hence for any k = 1, . . . ,m, 

N 



NlpkSf, A^^Sfc -pk 



(x + \J x + a 

by Theorem 1.1 in [5] when A^^ is given. Here the asymptotic variance can 
be computed using formula (4.23). From result (1) of Theorem 1.1 of [3] it 
can be concluded that 

Thus we are done by the Fubini theorem and the Cramer-Wold device. 



3. Proof of Theorem 1.2. Let z = u + iv, v > 0. Recall that the Stieltjes 
transform is defined for any distribution function F as 

mF{z)= f^—dF(x), zgC+ = {zGC,9z>0}. 
J X — z 

Hence, the Stieltjes transform of Gn{x) is 

1 ^ 1 

(3-1) '^GM{x)iz) = -^Yl 
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and it suffices to consider 

First we obtain a decomposition as follows: 



(3.2) 
where 



dHjx) 

bx — z 



K 



^1+1^2 + ^3 + ^4, 



1 



1 



1 ^ 

^£:tV(pfc/iV)trA^i-z (pfc/iV)trA-i 



1 ^ 



1 ^ 1 fdH{x) 



^ K — 2 J xb — z 

It is straightforward to verify that 



(3.3) 

then we have 



< 



1 



Pkb- z 



< 



K 



(3.4) 



E\VA< 



SfcA;- ^Sfc-^trA;,! 



N 



< 



M 



as iV — > 00. Similarly, by (3.3) one can find 



(3.5) 



E\V2\ < 



1 



K 



v^NK 



1 



v¥iK^^k 



< 



M 



k=i -+PfcS^A/sfc N 



as ^ 00 . From the uniform integrability of the random variable ^ tr A ^ , 
one can obtain 



(3.6) 



E\V-^\ < 



1 



K 



Y.PkE 



k=l 



-1 



It is obvious that IV4I converges to zero. This, together with (3.1)-(3.6), 
implies Theorem 1.2 and thus we are done. 
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4. Proof of Theorem 1.3. We begin the proof of this theorem with the 
replacement of the entries of S by truncated and centrahzed variables. Since 
Evfi < oo, we have e~^EviiI{\vii\ > e^/N) for any e > 0. Thus a posi- 
tive sequence £n converging to zero can be selected so that 

(4.1) e];,''Evtj{\vu\>eNVN)^0. 

Define Vij = %-/(|%-| < snVN) and Vij = Vij — Evij, i = 1, . . . ,N, j = 
1,. . . ,K. The corresponding matrices and vectors are denoted by 8^,8^, 
and Sfc, k = 1, . . . , K , the elements of which are Vij or Vij instead of Vij. 

Let 

A,, = SfcPfcS^ + a2l, Ak = SkPkSl + a^I. 
It follows from (4.1) that 

/ K _ K \ 

\fc=l fe=l / 

(4.2) <p(\j{\vij\>eN^)^ 



< NKP{\vii\ > EN^) 0, 



as A'" oo. 
Observe that 



(4.3) ^ 

" ^\\Eskf + 2\slAl^Esk\. 



a 

Concerning the first item on the right, we have 

(4.4) \\Eskf = \\Esif = o{N~'^ 
where we use the fact that 

(4.5) Evii = o{N^'^l^). 

For the second item, by (4.4) one can find 

1 



E{slA^^Esk? < —\\EskfE\\skf 



(4.6) 



1 



-^\\EsifE\\si\\^ = o{N-^'' 
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Combining (4.3)-(4.6), one can conclude that 



E 



(4.7) 



K 



K 



k 

Next we will show that 
(4.8) j:{slA^h,--slA-,h,)^0. 



By matrix inverse formula A 



-1 



B 



-1 



B^{B - A)A^\ we have 



ISfcA^^Sfc -s^A;, ^Sfcl 



(4.9) 



= \slA^^[{ESk)PkEiSl) 

- {ESk)PkSl - SkPk{ESl)]A^%\ 

<^\\skf\\ESkf\\Pk\\ 



+ \-siA^'[{ESk)PkSl + SkPkiESi)]A^'-Sk 



< 



MN. 



Sfcf (^{)n)^ 



+ \^^AI\E^^)VM + SfcPfcii;(SD]A^isfc| 
By Lemma 2.7 of [2] one can obtain 



si. A,T^SfcPfc(£;Si.)A,T^Sfc - -tr A^^SfcPfc(i?S^)A 



<^^£(i?T)n)^i?trS,S^ 



< 



[EvxxfEtY^x^l 



which implies 

K 



(4.10) ^ s^A^iSfcPfc(i?SDA^isfc - - tr A^iSfcPfc(£;sDA^i ^ 



fc=i 



Here we use (4.5) and 
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On the other hand we have 



ltrA^iSfcPfcii;(SDA,T^ 



\Evu 



Ar3/2 



< 



M\Evi 



K 



where e = (1, . . . , 1)'^. This, together with (4.5), gives 



E 



^EtrA^iSfcPfc(£;sDA^i 



k=l 



and by combining (4.10) one can then find 



K 



(4.12) 



^ slA^'S,Pk{ESl)A^'sk ^ 0. 



fc=i 



Similarly, one can also show that 

K 

(4.13) Y.^lA^\ESk)FkSlA^'s,^0. 

k=l 

Thus (4.8) immediately follows from (4.5), (4.9), (4.12) and (4.13). It is easy 
to check that 



(4.14) 



l-var(wii) = o(iV~^). 



Applying this and the argument similar to the centralization step, one can 
then renormalize the underlying random variables. Consequently, it can be 
assumed that the underlying random variables satisfy 



Ev 



11 



0, 



Ev 



11 



1, 



\vii\ < £nVN 



In the sequel we still use Vij, s^, Sk and A^ instead of vu, s^, Sfc and Afc 
to simplify the notation. 
Write Sk = s^A^^s^ and 

E Pk^k = E Pk^ki^k^ - A~^)sfe + tr A-^ PkSksl 

k=l k=l \ \k=l / / 



K 



V^^ + AT-a^trA-i 



k=i ^+PkSk 

Further, after some simple computations one can find 

K 



St 



1 



^^l+PkSk 



N-a^ tiA 



-1 



K. 
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Applying the formula 
(4.15) ^ 



PkSk - bNPk 



l+PkSk I + bNPk (1 +PfcSfc)(l +&ArPfc) 

to the above identity we can arrive at 



(4.16) 
where 



K 



(■X , U .2J2Pki^k-bN) = Ui-U2-U3 + Ui-U5, 



rj _^ Pl{Sk-bN? 



jj ^y- Pkisk-bNY 

^ ^^{'^+PkbN)^{l+PkSk)' 



Ui = N{1 - a^bN) -K + K J J 



and 



K 



U5 = "^Pkisk - Q 



k=l 



dHj\f{x) 
+ xbN 

dHiy{x) 



{l+PkbNr J i'^+xbN) 



As will be seen, the contributions from U3, U4 and C/5 can be ignored and 
the main terms are Ui and U2- 
It is easy to see that b]\f satisfies 



(4.17) 



1 

bN 



O + Cat 



xdH^ix) 
l + xb]\f ' 



that is, U4 = 0. 

For the term C/3 we have 



K 



m <MY,\sk- bNf < M(C/3i + C/32 + C/33) 

k=l 



where 



K 

k=l 
K 

C/32=E 
k=l 



Sfc- — trA^i 



1 *-i 1 * 1 
-trA,^--trA- 
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and 



U33 = K 



— trA-^-bN 



From Lemma 2.7 [2] one can find 



K 



k=l 

For the term C/32 



K 



k=l 



l+PkSk 



:0(Ar-2), 



where we use the fact 
(4.18) 



Pk^l^k^^k 



l+PkSk 



< 



In the sequel, we will not mention it again whenever (4.18) is used. By 
Theorem 1 of [3] we have 

^33^0. 

From the above argument it can be concluded that 
(4.19) 

We now analyze the term C/i by computing its variance: 



f4 2m f\ sr Pki^k-hN) 

^ ' ' \t, {1+bNPkr ti i'^+bNPk''^ 



^ piE{sk - h^y 

(1 



U11 + U12, 



where 



K 



V(l+ta.)3 {l+bNPkf 



and 



K 



Ui2= J2 E 

kij^k2 



vlXsk^-^N? pl^E{sk^-hNf 



{l + hNPk^T i^ + bNPki] 



p\{h2-bN? pl^E{sk2-hNf 



[l + hNPkiY iX + bNPkif 
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Similarly to the argument of (4.19), one can get 

K 

Uii<MY,E{h-hNf 

k=l 
K 



k=i 

4.21 

/I ^ ^ 

+ Mi^^; f —trA"^ -^AT 

= 0{eN), 

as N ^ oo. Indeed, in the last step we also use the fact that 

< (^l/a^ + bNy ^E{ti A"^ - Nbjvf ^ 0. 

To evaluate the term U12, we need to decompose it further as shown 
below. The strategy is to split A^^'^ into the sum of A^^^^^ and 

SO does for A^^^. Thus one can find 

(4.22) U12 = C/121 + f/122 + Ui23 + C/124 + • • • + Ui29, 

where 

K K 

Ul21= E E{pk^k2Pkik2Pk2ki), Ui22= E E{pkj^k2Ckik2Pk2ki), 

ki^k2 ki^k2 

K K 

Ui23= X! ^iPkik2Pkik2Ck2ki), ^124= ^ -^(pfcifca Cfcifc2Cfc2fci ) ) 

ki^k2 ki^k2 

K K 

^^125 = -2 E E{pk-i^k2Ckik20ik2ki), Ui2Q = -2 ^ -E^(j'fciA:2 /^fclfc2 "^2^1 ) ' 

ki^k2 ki^k2 

K K 

Ul27 = -2 ^ E{pk^k20tkik2Ck2ki)-, f^l28 = -2 ^ S(pfe^fc2 «fcife2/^fc2fcl ) ' 

fel^fc2 ki^k2 

K 

Ui29 = AY^ E{p 

kik2^kik2 ^k2k\ ) i 

kiy^k2 
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with 



Also, set 



Pij 



a 



[i + bNPi)Hi + bNPjy' 



E- 



(l+p,sjAris,)2 (l+p,sjAr.is,)2' 



7ii —PiP'j^i ^ij SjSj- A- - Sj, 



Pi7ii(sf A - .^Si - bj^) paij{sfA--^Si - bN) 



l+PjsjA^.hj 



E- 



l+PjsjA,/sj 



l+Pfc,s^^A^^Sfc/ 



"fc. 



1 +Pfc2Sfc2'^fclfc2^^2 

As will be seen, each of 3.nd Oij converges to zero in some way and the 

convergence rate is needed to attain our aim. In the subsequent paragraphs 
we show that each term Ui2j,j = 1, . . . , 9, converges to zero. 

Consider the term f/121 first. It is straightforward to verify that 

2 



E 



(4.23) 



T 1 

SfcBfcSfc - — trBfc 



1 ^ 9 

j^,{Evf,-3)Y^Eib^f + -^EtrB,B 



T 
k 5 



where = {b^jlj^) is any symmetric matrix independent of s^. It follows 
that 



^[(« A^/,^Sfc, - b^f - E{slAll,s,, - 6^)^)1 A^/,J 



and then, that 



U 



121 



E 



K 



X! Pklk2E{l3k^k2\-^klk2)^^l^k2kA-^klk2' 
■kxi^k2 



(4.24) 



E (i?(trA^.l-tri?A 



K 



k\k2 



2 x2 
kxk2' 



2\2x 
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Since the distribution of trA^^^^^ is dependent on different ki,k2, the 
difference between tr A^^^^,^ and tr A~^ caused by a different k must be ehm- 
inated. To this end, by sphtting A^^^^^ into the sum of A^^^ and ^fcifej, one 
can get 

(4.25) + ME{tr^k,k,A^^ - tvE^k,k,A^'f 

+ ME{trA^^-trEA^^f 

<M + M^(tr A^2 _ trEA^^f, 

where 

and we also used 

Repeating a step similar to (4.25), by Theorem 1 of [3] one can then conclude 
that 

kij^k2 

(4.26) < + __E(tr A-2 _ tri?A-2)2 

as ^ oo. Again, by an argument analogous to (4.25), one can find 

m^Klk2-^^^f 

< ME{tvA-^ - EtvA-^)"^ + M(^tr A"^ - NbN)^ + M. 

The second term on the right-hand side of the above inequality is bounded 
by the argument of Theorem 1 of [3]. We also claim that the first item 
on the right-hand side has an order 0{N). To see it, set J-j = cr{si, . . . , sj) 
and Ej{-) = E{-\J^j). By decomposing as the sum of a martingale difference 
sequence and using the Burkholder inequality, we have 

S(trA-^-^trA-^)^ 

(4.27) = E - Ek^i) tr(A-i - A^^) j 
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<ME 



<KMyE( (Ek-Ek-i) P^^k^k^^k 

ti ^ l + (pfc/iV)trA^i 



< KMY El sTEki ^ T Sfc 

fc=i ^ 



where 



^k — ^ , T A -1 ' 

and we also use (4.18) and the equaUty 

_ PfcSfc A^^Sfc _ VkiPkh - {Pk/N)tTA^.^) 
l + {pk/N)tvA^^ l + (pfc/iV)trA^i 

Combining the above one can conclude that 

^ E E{trA,^\^ - NbM)' < f - as iV - oo. 

ki=/=k2 

The basic inequality < EX^ implies the remaining term in (4.24) 

also goes to as ^ oo. Hence f/121 can be ignored. 
Consider the term U122 second. By (4.23) we have 

E{{pk2 afc2 Sfc^ A^ ^ Sfc2 Sfc2 A^ ^ s^, ) | s^^ , A^ ^ ) 
(4.28) ^ 
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It follows that 

K 

Ui22= 5Z -^[-^bfcifc2Cfcifc2|Sfc2; A^ifcJAafci] 
M ^ 



^ M ^ ^( \ .1 , \ 



k\i=ki 

MK^ MK^ , 1 ,2 



Similarly, the term U123 converges to zero. 

Third, consider the term ^7i24- To simplify the notation, we write 

^ki — ^ki-^k^k2^'^i ' ^k2 ~ ^k2-^kik2^k2 ' J = 1) 2, 

Skik2 = ^ki-^kik2^k2^k2-^kik2^^'i ■ 

According to (4.28) one can find that 

K 

Ul2A= J2 EiPkik2(yi^lkik2lk^k2(^\) 
k\i=-k2 

y 2 2 771/ (2)n2c./ (2)x2 

(4.30) 

<M ^ Eal^{pk2Sk,k2) + 

k\i^k2 

<MeN = o{l). 



J 



Fourth, since the composition of the terms U125 and Ui27 is similar, we 
analyze only the U125 term. From (4.28) we obtain 

K 

(4.31) C/i25 = -2 ^ -E^(Pfcifc2«i2"fci7feifc2Pfc2 7fcifc2(Sfc2^ -&iv)) 



A' 



('^•32) +772 E l'fclfc2^(PfciPfc24?"fc2)^^(Pfc2 7fcifc2(4? - 

fei^fc2 
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By the Holder inequality and (4.28) one can find 

l^(Pfc2 7fcife2(Sfe - &7v)afci)| 
(4.33) < M{E{pu,h,k,o.u,ffl\E{s'^^} - hM^f'^ 



< 



M 



which implies that the term in (4.31) converges to zero. Moreover, the ab- 
solute value of each summand in (4.31) is not larger than 



ME\al {pkjk^kifis'j}^ - bN)\ 



ME[\4^ - bN\E{\pkJk^k,faljA^^^^,Sk,)] 



<ME 



Ji) 



bN\E 



1 (2) 



«fc2 



M 
iV3 



E[\s'iybME{\pkJSfal\A^^,^,Sk, 



Ar3 



< 



M 



which leads the sum in (4.31) to converge to zero. So C/125 converges to zero 
as ^ 00. 

Fifth, consider Uue (C/128 can be analyzed similarly): 



26 



K 



(4.34) 



(4.35) 



-2 Pkik2E[pl,{s)}^ -bNfpk2lkik2{Sk^ -bN)aki] 

kij^k2 



K 



+ 2 51 Pkik2PkiE{s[^' - In) Epk^-ik^k2{sl2 -^N)aki- 

kiy^k2 



From the estimate (4.33), the sum in (4.35) has an order 0{N-'^^'^). On the 
other hand, each summand in the sum in (4.34) can be rewritten as 

E[p\pk2Pk^k2{Sk} - bNfak^E{-fk^k2{Sk2 - bN)\A.^X^,Sk^)] 
<E[plakAs^^ -bN)\E{^l,JAl^\^,Sk,)f' 
x(E((.«-6^)V.-\,,s,J)V2] 
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M 
- Ar5/2' 

and, thus, the sum in (4.34) converges to zero. Thus U126 converges to zero 
as well. 

Finally, consider the term C/i29- As was done for the terms C/124, U125 and 
Ui2Q, the C/129 term is split into the sum of two terms. It follows from (4.33) 
that one of them, 

K 

fcl^fe2 

converges to zero. The other term is 

K 

4 H ^(Pfclfc2^'fcl7fclfc2(4? -^A^)«fc2Pfc2 7fcifc2(Sfe -&iv)«fcl)• 
The absolute value of each of the above summands is not larger than 

{E{pk^Skik2{Sk} -hN)ak^ff''^ X {E{pk^Sk^k2{s^k2 - bN)ak2f f'^- 
Note that 

E{pkiSkik2{Sk} - bN)akjf 

= E{{{s'i^ - hN)ak,fE{{pkJk,k2f\Klk2^^k,)) 

= ]^^((4? - ^^)pfci4?«fci)^ 

Similarly, 

E{pk2Sk,k2{4'2 - bN)ak2f < M/N\ 

Hence C/129 converges to zero. 

Summarizing the above argument, one can conclude that the variance 
of the term Ui converges to zero as iV — > 00 and thus, it is sufficient to 
compute the asymptotic value of its expectation, which can be accomplished 
as follows: 

K 

YpkE{sk - hNf 

k=l 



fc=l 

K 
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(4.36) =Y.PkE(h -^tixA^ + Y.p^e( -^trAl^ - ' 

k=i \ / 

2 ^ 

= l|2^trA-2^Pfc + o(l) 

A:=l 

^2c / y / -:rdH{x), 

J (x-h 0-2)2 y (l-^2;6)3 ^ ^' 

where =p|/(l -|- bj\fpk)^, and in the second and third equahties we use a 
trick similar to (4.25). 
Note 

(4.37) F'^^^ = il-c)I[o,oo) + cr''', 

where F^'^'^ represents the hmiting spectral distribution of P"'^/2s^spi/2 
with P^/2 = diag(y/pi, . . . , y/px). From (4.37), one can get 



tr - NbN = tr(pi/2S^Spi/2 + a^iy^ -K j 



dF^N,HN(^X^ 

X + a"^ 



and Theorem 1.1 of [3] is then applicable. Thus, for U2 one has a central limit 
theorem and it then suffices to show that converges to zero in probability. 
It is obvious that the term C/5 becomes zero when pi = • • ■ = pk] however, its 
convergence in probability appears to be somewhat troublesome when the 
powers of the users are not the same. We will provide an abridged analysis 
for this case. Set 



'^fc — 7T~; — ; — v7> ^ 



dHj\[{x 



(1+Pfc6^)2' J (l + x6^)2- 

Using steps analogous to (4.25), one can obtain 

K 

I ,S;, — 
.=1 ^ ^ 

K 



+ (trA ^ -NhN)^ Y.Pk{ak-a) + Op{l) 



^k=i 



^Pk(h- j^'^T:A^^]{ak- a) + Op{l). 



k=i 
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Hence it is necessary to show that 



I K , , 
(4.38) iJ^tE\^pk\sk-^^^K^ 



(flfc - a) 



converges to zero. Expanding out the right-hand side of (4.38) one can get 



where 



t^5i = X^£^(^Pfc(^Sfc - -^tr A^.^^(afc - a)^ , 



ki^k2 



It is easy to see that 



X I s^k2^ - ^ tr A {ak, - a^ak^ - aU . 



M 2 
U51 < — 2^{Pk{ak - a)) ^0. 

k=i 

Regarding the term Ufj2, one can show that it converges to zero by an argu- 
ment similar to that used for the preceding term U12 and since the process 
is somewhat tedious, it is omitted. 

For the computation of (1.11) and (1.12), without loss of generality, sup- 
pose p = I; otherwise replace cr^ by cr^/p. As for the formulas (1.9) and 
(1.11) one can refer to, respectively, (1.18) and (5.13) of [3]. 

Now let us derive (1.12). It is shown in [7] that m(z) = mpc,H{z), for each 
z £ C"*", is the unique solution in C"*^ to the equation 

From this equation, the inverse function has an explicit form 

1 ftdHlt) 

(4.40) z = + c ' ^ ' 



m J 1 + tm 
and one can then find for H{t) = /[i,oo)(i) 

(4.41) z = ^ + 



m{z) 1 + m{z) 
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Suppose the m2 contour encloses the mi contour (see [8] or [3] for the range 
of m(x) for a real x and contour of m). For a fixed m2 it follows from (4.41) 
and the Cauchy residue theorem that 

dmi 1 f {ml + mi) drill 



r ami - ^ [ 

J izimi) + a'^){mi — m2Y cr"^ J 



{z{mi) + a'^){mi — m2y cr'^ J {mi — ma){mi — mb){mi — m2y 



2ni m? + mi, 



(m6-m2)2(m6-m„)' 

where 



Consequently, 



_ -(1 + (c - + v/(l + (c - 1)/(T^)^ + 4/^ 

ma- 2 

-(1 + (c - 1)/ct2) - V(l + (c-l)/a^)^ + 4/a^ 



-I- mb 



^ a'^ni J {m^ — m2Y{mb — ma){z{m2) + o'^) ^ 
m2+m, /• (m2 + m2)dm2 



a^TT{mb-ma)i J {m2 - m^)^ {m2 - ma) 
2{ml +mf,){ml + ma) _ 2c 



a^{mb-ma)^ ((a^ + c - 1)2 + 4cj2)2 ' 

5. Proof of Corollary 1.1. Using the Taylor expansion, one can find 

K 

Y,{l0g{l + f3k)-\0g{l + bNPk)) 
k=l 

^ Pk - bNPk _ y> {Pk - ^NVk)^ y> {Pk - bNPk)^ 
t^l^ + bNPk tl'^{l + hNPk)^ fctl 3(1+6)^ 

dHN{x) ^ ^ {f3k - bNPk)'^ , ^{Pk-bNPk) 



f dHN[x) ^ M-^(Pk-bNPk) ^Y^lPfc- 

-|:(A--)(TTk-/^)' 

where each is located in the interval [Pk , ^NPk] ■ Since 
^ (/3fc - bNPk)^ ^^\a L |3 

Corollary 1.1 holds by the argument of Theorem 1.3. 
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